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ABSTRACT 

We report on a multiwavelength campaign on the radio-loud Narrow-Line Seyfert 1 (NLS1) 
Galaxy PMN J0948+0022 (z = 0.5846) performed in 2010 July-September and triggered by 
a high-energy 7-ray outburst observed by the Large Area Telescope (LAT) onboard the Fermi 
Gamma-ray Space Telescope. The peak flux in the 0.1 — 100 GeV energy band exceeded, for 
the first time in this type of source, the value of ~ 10~ 6 ph cm~ 2 s _1 , corresponding to an 
observed luminosity of ~ 10 48 erg s _1 . Although the source was too close to the Sun position 
to organize a densely sampled follow-up, it was possible to gather some multiwavelength data 
that confirmed the state of high activity across the sampled electromagnetic spectrum. The 
comparison of the spectral energy distribution of the NLS1 PMN J0948+0022 with that of a 
typical blazar - like 3C 273 - shows that the power emitted at 7 rays is extreme. 

Key words: galaxies: individual: PMN J0948+0022 - galaxies: jets - galaxies: Seyfert 



1 INTRODUCTION 

Relativistic jets are the most extreme expression of the power 
that can be generated by supermassive black holes at the centres 
of galaxies. Their bolometric luminosity can reach values up to 
~ 10 49-50 erg s _1 , peaking in the 7-ray (MeV-TeV) energy band. 
To date, two classes of active galactic nuclei (AGNs) are known to 
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generate these structures: blazars and radio galaxies, both hosted in 
giant elliptical galaxies (Blandford & Rees 1978). The first class is 
composed of AGNs with the jet viewed at small angles and, there- 
fore, the effects of the special relativity play a dominant role in 
shaping the emission of the electromagnetic radiation. In the sec- 
ond class - radio galaxies - the jets are viewed at larger angles and 
therefore the Doppler boosting is less intense. 

The recent discovery of variable 7-ray emission from 4 NLS Is 
confirmed the presence of a third class of AGNs with relativistic 
jets (Abdo et al. 2009a,b,c, Foschini et al. 2010a), which was al- 
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Figure 1. Fermi/LAT lightcurve (0. 1-100 GeV) of PMN J0948+0022 in 2010, binned over 1 month during the first six months of 2010 and with 3 days time 
bin for the period 2010 July-September. The value on the abscissa corresponds to 2010 June 30 (MJD 55377). Blue points are detection with Test Statistic 
TS > 10 (equivalent to ~ 3cr, see Mattox et al. 1996 for the definition), while the others are 2cr upper limits. The abscissa has the same scale as Fig. [2] to 
make comparison easier. 
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Figure 2. Lightcurve of PMN J0948+0022 at optical (ATOM, R filter, top panel) and radio frequency (OVRO, 15 GHz, bottom panel). The abscissa has the 
same scale as Fig.[T] to make comparison easier. The value corresponds to 2010 June 30 (MJD 55377). 
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ready suggested by the detection of variable radio emission with flat 
or inverted spectrum (e.g. Grupe et al. 2000, Komossa et al. 2006, 
Doi et al. 2006, Yuan et al. 2008). This finding poses intriguing 
questions about jet systems and how these structures are generated. 
One of these sources, PMN J0948+0022 (z = 0.5846) is classified 
as a typical NLS1, i.e. FWHM(H/3)< 2000 km s~\ [OIII]/H/3 < 
3, with the presence of the Fell bump (see Pogge 2000 for a re- 
view on the characteristics of NLSls). However, it displays also 
strong, compact and variable radio emission, with inverted spec- 
trum, suggesting the presence of a relativistic jet (Zhou et al. 2003). 
The confirmation came with the detection of high-energy variable 7 
rays by FermifLAT (Abdo et al. 2009b, hereafter "Discovery 2008"; 
Foschini et al. 2010a). A multiwavelength (MW) campaign per- 
formed in 2009 March- July displayed coordinated variability at all 
frequencies, thus confirming that the source detected by FermifLAT 
is indeed the high-energy counterpart of PMN J0948+0022 (Abdo 
et al. 2009c, hereafter "MW Campaign 2009"). However, the vari- 
ability observed in 2009 was much smaller in amplitude (more than 
a factor 2 at 7 rays) with respect to the 2010 campaign. 

The main differences between NLSls and the other AGNs 
with jets are the optical spectrum and the radio morphology, which 
is quite compact and without extended structures (Doi et al. 2006, 
Komossa et al. 2006, Abdo et al. 2009c, Foschini et al. 2010a, Gu 
& Chen 2010). In addition, NLSls are generally hosted in spiral 
galaxies (Deo et al. 2006, Zhou et al. 2006), which is at odds with 
the current observational paradigm linking very powerful jets with 
ellipticals (cf, e.g., Marscher 2009; see also Bottcher & Dermer 
2002). Obviously, to make a proper comparison we would need ei- 
ther a direct imaging of the host galaxies of the four 7-NLSls or 
the morphologies of a similar sample in term of redshift and mag- 
nitude distribution. Our sources have redshifts between 0.061 and 

0. 585, while the work by Deo et al. (2006) is dedicated to a sample 
with 2 < 0.03 and Zhou et al. (2006) studies the morphology of 
NLSls with z < 0.1, so there is only a partial overlap. At least one 
of the four NLSls detected at 7 rays is definitely hosted by a spiral 
galaxy (1H 0323+342, see Zhou et al. 2007, Anton et al. 2008), 
while for the others there are no high-resolution observations avail- 
able. Therefore, some halo of doubt remains. 

Despite these uncertainties - worth studying - there are other 
novelties introduced by the 7-ray detection of NLSls. The spectral 
energy distributions (SEDs) of the four 7-NLSls point to AGNs 
with masses of the central black hole in the range of 1O 6 - 8 M 
and high accretion rates (up to 90% of the Eddington value). These 
values are common for NLSls (e.g. Boiler et al. 1996, Grupe 2004, 
Collin & Kawaguchi 2004), but not for blazars or radio galaxies. 

Another piece of the puzzle is given by the search for the par- 
ent population. Indeed, if 7-NLSls seem to be similar to blazars, 

1. e. with the relativistic jet viewed at small angles, there should be 
a parent population with the jet viewed at large angles (as in the 
case of blazars vs radio galaxies). The first source of this type - 
PKS 0558-504 (z = 0.137) - has been recently found by Gliozzi 
et al. (2010). Therefore, it seems that NLSls could be a set of low 
mass systems "parallel" to blazars and radio galaxies. 

A key unknown in this research field was the power re- 
leased by jets of 7-NLSls. During the early observations and 
the 2009 MW campaign the maximum measured luminosity of 
PMN J0948+0022, the most powerful of these NLSls, was ~ 
10 47 erg s _1 (0.1-100 GeV). On the other hand, it is known that 
blazars can reach greater luminosities (see, for example, Ghisellini 
et al. 2010), up to ~ 10 49 erg s _1 in the case of the bright- 
est outbursts of 3C 454.3 (Bonnoli et al. 2011, Foschini et al. 
2010b) or even in excess of 10 50 erg s _1 in the recent episode 



Table 1. Summary of the spectral fitting of the Fermi/LAT data integrated 
on long timescales. -Fe>ioomcV is in units of [10 — 6 ph cm -2 s -1 ]. The 
statistical measurement errors are at la. Systematic errors are not included. 
The latest estimates for the adopted Instrument Response Function (IRF) 
P6.V3.DIFFUSE are: 10% at 100 MeV, 5% at 500 MeV and 20% at 10 
GeV, while the error in the photon index is ~ 0.1 (Rando et al. 2009). 
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0.23±0.01 
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Jul 07 - 10 
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Figure 3. High-energy 7-ray SED of PMN J0948+0022. The red star sym- 
bols refer to the period of maximum flux (2010 July 7—10, 4 days), the 
blue squares indicate the emission integrated over the period 2010 Au- 
gust 1 — September 15 (45 days), while the magenta circles are for the 
period 2010 June 1—30 (30 days). The grey triangles refer to the quies- 
cence period of one month centered on 2008 December 5 (corresponding to 
the Swift observation, Abdo et al. 2009b). All the detections have TS > 10 
(~ 3<t, see Mattox et al. 1996). Upper limits are at 2cr level. See also Ta- 
ble 1. 



of the gravitationally lensed quasar PKS 1830—21 (Ciprini et al. 
2010). The question of the power could be studied in 2010 July, 
when PMN J0948+0022 underwent a period of intense activ- 
ity (Donato et al. 2010, Foschini 2010) with a peak flux of ~ 
10 -6 ph cm -2 s _1 (0.1 — 100 GeV), corresponding to an observed 
luminosity of ~ 10 48 erg s . Here we report about the study of 
this period. Some very preliminary results have been presented in 
Foschini et al. (2010c). 

In the following, we adopted a ACDM cosmology as calcu- 
lated in Komatsu et al. (2009) based on WMAP results: h = 0.71, 
Sl m = 0.27, Ha = 0.73 and with the Hubble-Lemaitre constant 
H = lOO/i = 71 km s _1 Mpc -1 . 



2 DATA ANALYSIS 

The analysis of the data from the different instruments and facili- 
ties participating in the present campaign was done in basically the 
same way as for that of 2009 (see MW Campaign 2009 paper and 
references therein). There were some updates in the software used: 
SwiftfXKT and UVOT data were studied with the specific packages 



4 L. Foschini et al. 



Table 2. Summary of results from analysis of the Swift data obtained on 2010 July 3 (ObsID 00038394002) compared with the observation of 2008 December 5 
(ObsID 00031306001, Discovery 2008) and 2009 May 15 (ObsID 00031306006, MW Campaign 2009), both reanalyzed here with the latest software (with 
consistent results), when the source was at low 7-ray flux. The value of Njj was fixed to the Galactic column as measured by Kalberla et al. (2005). 



XRT (0.3 - 10 keV) 


ObsID 


Date 


Exp. 


N„ 


r 


Flux .3_l0kcV 


X 2 /dof 




[YYYY/MM/DD] 


[ks] 


[10 20 cm" 2 ] 


[10 


~ 12 erg cm -2 s _1 ] 




00038394002 


2010/07/03 


1.6 


5.22 


1.5±0.3 


5.5±0.5 


0.53/2 


00038306006 


2009/05/15 


1.3 


5.22 


1.7±0.3 


2.3±0.3 


(*) 


00031306001 


2008/12/05 


4.2 


5.22 


1.8±0.2 


3.9±0.2 


3.2/8 


Low statistics. The 


fit has been performed by u: 


;ing the maximum 


likelihood (Cash 1979), obtaining c-stat of 57 for 58 d 


.o.f. 




UVOT (observed magnitudes) 


ObsID 


V 


6 


u 


uvwl 


uvm2 


uvw2 


00038394002 


17.6±0.2 


18.3±0.1 


17.24±0.08 


17.19±0.07 


17.34±0.09 


17.23±0.06 


00031306006 


18.3±0.2 


18.4±0.1 


17.7±0.1 


17.43±0.08 


17.5±0.1 


17.65±0.07 


00031306001 


18.2±0.1 


18.56±0.08 


; 17.79±0.07 


17.57±0.06 


17.63±0.07 


17.63±0.05 



included in HEASOFT V. 6.1C0 and the calibration database up- 
dated on 2010 September 30. 

The light curve obtained from FermifLKT data is shown in 
Fig.Q] while the spectra are displayed in Fig.|3]and the fits are sum- 
marized in Table 1 . The change in the photon index, confirmed by 
the spectra shown in Fig. [5] indicates some "harder when brighter" 
behavior. This comparative analysis should be considered with care 
and the caveats about the comparison of fits over the whole energy 
band versus fits on separate energy bins are described in detail in 
the Section 4.4 of Abdo et al. (2010a). 

Fig. E] shows the optical (R filter) and radio (15 GHz) light 
curves from ATOM (Automatic Telescope for Optical Monitoring 
for H.E.S.S., Namibia) and OVRO (Owens Valley Radio Observa- 
tory) 40 m telescope, the latter as part of an ongoing Fermi blazar 
monitoring program (Richards et al. 2009). Multifrequency obser- 
vations to measure spectral changes were also done at Effelsberg 
(Project F-GAMMA, Fuhrmann et al. 2007) and RATAN-600 (Ko- 
valev et al. 1999), respectively. The study of the morphology was 
done with Very Long Baseline Interferometry (VLBI) observations 
performed on 2010 September 17 as part of the MOJAVE mon- 
itoring prograrr0 conducted with the Very Long Baseline Array 
(VLBA) at A = 2 cm and provided high resolution total inten- 
sity and linear polarization images. Detailed description of the MO- 
JAVE program as well as method of observations, data processing 
and analysis is presented by Lister et al. (2009). 

The source was close to the Sun in 2010 July and there- 
fore it was not possible to perform a densely-sampled optical-to- 
X-ray follow-up. Only one X-ray observation close to the out- 
burst remained and was performed by Swift on July 3 (Ob- 
sID 0003839400©, about 4 days before the outburst. A summary 
of XRT and UVOT analysis is displayed in Table 1, compared with 
the observation of 2008 December 5 (ObsID 00031306001, Dis- 
covery 2008) and of 2009 May 15 (ObsID 00031306006, the low- 



1 Including the XRT Data Analysis Software (XRTDAS) developed under 
the responsibility of the ASI Science Data Center (ASDC), Italy. 

2 http : / /www .physics . pur due . edu/astro/MOJAVE/ 

3 Within the project Swift-XRT Monitoring of Fermi-LAT Sources of In ter- 
est. See: http : / / www . swift. psu. edu /monitoring/. 



est X-ray flux measured during the MW Campaign 2009), when the 
source was in a low activity state at 7 rays. 

The data available confirm the outburst first detected at 7 
rays of the NLS1 PMN J0948+0022. The optical light curve from 
ATOM (R filter) indicates an increase of the optical activity with a 
change of about one magnitude during the early six months of 2010 
(Fig. [2] top panel). This trend is confirmed by the Swift/WOT ob- 
servation performed on July 3, when compared with an observation 
of the source at low 7-ray state (2008 December 5 in the Discov- 
ery 2008; 2009 May 15 in the MW Campaign 2009): the decrease 
of the observed magnitude ranges from 0.6—0.7 in the v filter to 
0.3—0.4 magnitudes in the ultraviolet filters. The same behavior 
was observed during the MW Campaign of 2009 March- July, when 
the source flux was decreasing, but showing stronger variability at 
optical frequencies rather than at the ultraviolet ones. This is likely 
to be due to the fact that ultraviolet frequencies sample the emis- 
sion from the accretion disk, while the optical emission is mainly 
due to the synchrotron radiation. 

The X-ray flux as measured by Swift/XKI was ~ 40% greater 
than the value measured in the Discovery 2008 paper and a factor 
~ 2.4 greater with respect to the lowest flux measured during the 
MW Campaign 2009. The low statistics in the short observations of 
2010 July 3 and 2009 May 15 prevent any conclusion on possible 
spectral changes. 

At 7 rays there was the most dramatic change in the emission 
of the electromagnetic radiation with an increase of flux by a factor 
~ 5 in the 0.1 — 100 GeV energy band. The outburst (Fig.[T} was 
characterized by a first sharp increase of the flux on July 7 (Donato 
et al. 2010) with e— folding timescal^fl of 3.6± 1.8 days, reaching 
the peak value of ~ 10 -6 ph cm -2 s _1 and remaining at such level 
for about four days. On July 1 1 the source was not detected (with a 
2cr upper limit of 2.7 x 10~ r ph cm~ 2 s" 1 ) and then it was detected 
again on July 12 above ~ 10 -6 ph cm -2 s _1 (e— folding rise time 
< 1 day). After a few days of low activity, there was a second 
shorter outburst (Foschini 2010) with flux ~ 10 -6 ph cm -2 s _1 



4 Some information were derived from the analysis of the LAT light curve 
with 1 day time bin not shown here. 

5 That is the characteristic time r to have an exponential increase or de- 
crease of flux: F(t) = F(t ) exp[-(i - t )/r\. 
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Figure 4. Total intensity and linear polarization images observed by VLBA 
at 15 GHz in different epochs as part of the MOJAVE program. Naturally- 
weighted total intensity images are shown by black contours, the con- 
tours are in successive powers of two times the base contour level of 
0.2 mjy beam" 1 . Electric polarization vectors direction is indicated on 
the left hand side by blue sticks, their length is proportional to the polar- 
ized intensity. Linear fractional polarization is shown on the right hand side 
overlaid according to the color wedge. 
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Figure 5. Spectral evolution at radio frequencies as measured by RATAN- 
600 (2010 August 13—29) and Effelsberg (all the remaining measure- 
ments). 



and rise time < 3 days. Then the source slowly returned to quies- 
cence, with some occasional detection at moderate flux. 

The spectral behavior is more difficult to study, as previously 
outlined. Nevertheless, the comparative analysis based on long 
timescale (month) integration, thus with a more robust statistics 
than that available on short timescales (Fig. [3] Table 1), indicates 
that there is some spectral hardening between the spectrum mea- 
sured during the four days of the peak and the spectra measured on 
the data integrated on 2010 June and August-September. It is not 
possible to firmly establish if this is due to a change of the peak of 
the inverse-Compton emission, as occurred in the case of 3C 454.3 
(see Bonnoli et al. 201 1), or if it is due to an effective increase of 
the population of GeV photons. 

The spectral evolution of the radio emission, as measured 
by Effelsberg and RATAN-600 (Fig. [5} indicate instead a pro- 
gressive inversion of the spectral index between 4.8 and 10 GHz 
(S v oc i>- a ), changing from a = 0.02 ± 0.03 on 2010 May 1 
(MJD 55317.83) to a = -0.81 ± 0.02 on 2010 June 26 
(MJD 55373.77), an average a = -0.80 ± 0.21 measured in the 
period 2010 August 13—29, and up to a value of a — —1.0 ± 0.1 
measured on 2010 September 18 (MJD 55457.45). 

We compare these results with the MW Campaign of 

2009 March-July, reported in Abdo et al. (2009c). In that case, 
the source was at a moderate high flux in 2009 April, peaking on 
April 10 with a flux of (4.2 ± 0.9) x 10" 7 ph cm" 2 s" 1 in the 
0.1-100 GeV energy band and T = 2.1 ± 0.1 (TS = 46; rean- 
alyzed in the present work). However, the flux integrated over the 
whole month of April gave a softer spectrum (r = 2.7 ± 0.2, cf 
Abdo et al. 2009c). In the 2010 July outburst, the 7-ray spectrum 
seems softer before (June) and after (August-September) the burst, 
with some hardening during the peak of the emission (July). The 

2010 July burst was a factor of ~ 2 greater than the 2009 April 
flare and also longer (more than 4 days vs 1 day). 

It is also interesting to note the evolution of the radio emission, 
polarization and morphology as shown by the VLBA observations 
at 15 GHz (Fig. |4). The observation performed on 2010 Septem- 
ber 17 shows a bright unresolved core with total parsec-scale flux 



° As indicated in the report #44 (2009 April 1 —7) of the Fermi Gamma ray 
Sky blog http : / / f ermisky .blogspot . com/. 
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Figure 6. (left panel) SED of PMN J0948+0022. The grey symbols refer to the low activity period observed by Swift (5 December 2008) and Fermi (one 
month centered on the Swift observation). The red symbols, used to model the SED, refer to July 7—10 (Femw'/LAT), July 3 (Swift XRT and UVOT), and 
June 26 (ATOM R filter). The model components are: synchrotron (green dashed line), accretion disk and torus (black dotted line), synchrotron self-Compton 
(grey long-dashed line), external Compton (black dot-dashed line) and the sum of all of them (blue continuous line), (right panel) SED with radio data from 
RATAN-600 (August 13—29, green open triangles) and Effelsberg (September 18, black filled circles). The orange continuous lines are the result of the 
modeling of radio data (synchrotron plus external Compton with seed photons from the molecular torus). The dashed orange line is the emission from the 
torus and the accretion disk. The blue continuous line and the red points are the SED of the peak (July 7—10) shown in the left panel. 



density as high as 538 mJy and a high degree of polarization 
(3.3%), greater than the values measured in 2009 (< 1%). Specifi- 
cally, when studying the evolution from 2009 May (MW Campaign 
2009) to 2010 September, we noted a drop in the polarized flux den- 
sity in the third epoch, coupled with a swing of the electric vector 
polarization angle (EVPA) of about 90°, which is maintained in the 
last epoch, even with a much greater polarized flux density (a factor 
~ 6). A similar behavior has been observed - for example - before 
the 2008 August outburst at 7 rays of the flat-spectrum radio quasar 
PKS 1502+106 (z = 1.839, Abdo et al. 2010b). 

It is not possible to estimate the speed of the radio jet on the 
basis of these 4 epochs only, particularly as in the present case when 
the source is strikingly compact. Observations are continuing, and 
we anticipate having a robust speed estimate in the near future using 
a larger sample of epochs. 



3 MODELING OF THE SED 

The spectral energy distribution shown in Fig. [6] (left panel) 
has been fitted with the model described in Ghisellini & Tavec- 
chio (2009), the same used in the previous works (Abdo et al. 
2009a,b,c). The electromagnetic emission from a relativistic jet 
is modeled with different components: the synchrotron radiation 
plus the inverse Compton emission derived by the interaction 
of the relativistic electrons with different types of seed photons 
(synchrotron, molecular torus, broad-line region, accretion disk, 
corona). We present here our favored interpretation of the collected 
data. Changes in the model parameters are negligible with respect 
to variations of ±30% in the 0.1-100 GeV flux or ±0.1 in the 
photon index (see Section 4 of Discovery 2008 paper for a more 



extended discussion on the impact of measurement errors on the 
calculations of the model parameters). 

The 7-ray flux is most important, because it flags the jet power. 
The extreme Compton dominance suggests a Doppler factor greater 
than that measured during the MW Campaign 2009: therefore, we 
adopted in the model a bulk Lorentz factor V = 16 (vs 10 of the 
MW Campaign 2009), with a jet viewing angle of 3° (vs 6°). This 
is one "economical" possibility, in terms of energetic. It would be 
possible to obtain the same Compton dominance with a lower mag- 
netic field, but with greater injected power (see below). In addition, 
it is worth noting that the magnetic field is linked to the X-ray emis- 
sion through the synchrotron self-Compton process (it is peaking at 
these energies) and the only available X-ray measurement has been 
performed well before the beginning of the outburst (Table 2). 

The energy distribution of the injected electrons (with power 

2.1 x 10 43 erg s _1 in the comoving frame) spans from 7 c , m in = 1 
to 7c,max = 2700 and is described by a broken power-law model 
with index -y% and 7,T 3 ' 2 below and above the break 7 C . break = 
200, respectively, where 7 C indicates the random Lorentz factor of 
the electrons. The external Compton peak is produced by electrons 
with 7c, P cak — 170, scattering off seed photons of the broad-line 
region. 

Our model shows that most of the dissipation occurs at about 

1.2 x 10 17 cm (2600 Schwarzschild radii, within the broad-line 
region), the magnetic field is B — 2 G. The mass of the cen- 
tral black hole was fixed to the value of 1.5 x 10 s A/0, as mea- 
sured by fitting of the accretion disk emission (see Discovery 2008) 
and in agreement with the measurement by Zhou et al. (2003) 
by using FWHM of emission lines. The accretion rate is at 50% 
of the Eddington value. The calculated jet powers (the output of 
the model) are Pradiat™ = 5 x 10 45 , P ma gnetic = 1-2 x 10 45 , 



First j-ray outburst of the NLS1 PMN J0948+0022 7 



^electrons = 6 x 10 44 and Pprotons = 10 47 erg s _1 , by assuming 
one proton per electron. 

These values have to be compared with those reported in pre- 
vious studies (specifically, see Table 4 of MW Campaign 2009 pa- 
per). The present values of the jet powers are generally greater, 
up to one order of magnitude in the case of the protons. Other 
differences with respect to the previous works are a lower mag- 
netic field with respect to the average of the MW Campaign 2009 
(B = 4.1 G), but similar to the value in Discovery 2008 (B — 
2.4 G). The injected power is similar to the previous studies, but the 
shape is different, with a greater 7 e ,max (2700 vs 2000 and 1600, 
in the MW Campaign 2009 and Discovery 2008, respectively) and 
a smaller 7 c , b rcak (200 vs 530 and 800, in the MW Campaign 2009 
and Discovery 2008, respectively). 

Since optical-to-X-ray data have been measured a few days 
before the outburst, it is possible that they are underestimated. If 
we assume an increase by a factor 2-4 in the optical-to-X-ray fluxes, 
keeping the 7-ray flux to the observed peak value, then it is possi- 
ble to model the SED with a relatively smaller bulk Lorentz factor 
(T = 13). The smaller beaming is compensated by an increase of 
the injected power of about one third, plus a greater magnetic field 
(B = 3.6 G). The jet powers (protons, electrons, ...) have negli- 
gible changes, because they are strongly dependent on the 7-ray 
flux. 

The region emitting the optical-to-7-ray flux is too compact 
to be responsible for the radio emission (its synchrotron flux is 
self-absorbed). This is confirmed by the variability analysis: the 
timescales at 7 rays are short, of the order of days or even less than 
one day, thus requiring a very compact source. About two months 
after the outburst at 7 rays, the radio flux density reached its max- 
imum (Fig. [2] bottom panel). During this period, the emitting re- 
gion can move outward and expand, thus becoming optically thin 
at radio frequencies. Moreover, near-simultaneous OVRO and MO- 
JAVE observations indicate < 50 mJy of arcsecond-scale flux den- 
sity at 15 GHz, hence nearly all of the radio emission in the source 
is generated on pc-scales. Therefore, to account for the radio emis- 
sion, we considered an additional larger emitting zone (Fig. [6] right 
panel). For the sake of simplicity, we assume the same bulk Lorentz 
factor of the optical-to-7 rays fit (V = 16). This region is ~ 0.6 pc 
in size, about two orders of magnitude larger than the 7-ray emit- 
ting region, and farther (5.8 pc) from the central black hole. At this 
distance, the seed photons for the external Compton process are 
from the molecular torus, emitting at infrared wavelengths. 



4 DISCUSSION AND CONCLUSION 

We presented the results of a MW Campaign from radio to 7 rays 
performed in 2010 to study the evolution of the electromagnetic 
emission of the Narrow-Line Seyfert 1 Galaxy PMN J0948+0022. 

Despite all the caveats and uncertainties, what clearly emerges 
from this episode is the ability of a NLS1 to generate an extreme 
power (~ 10 4S erg s _1 ), when compared to that of a typical flat- 
spectrum radio quasar. Fig. |7J shows a comparison of the SED 
(rest frame) of PMN J0948+0022 with that of the archetypical 
blazar 3C 273, known to have a strong disk emission too. At first 
look, the figure shows that the NLS1 has a more extreme Comp- 
ton dominance: although its radio-to-X-ray luminosity is smaller 
than that of 3C 273, the 7-ray power is greater. If we renormal- 
ize the two SEDs by taking as reference the peak due to the ac- 
cretion disk emission, then it is necessary to multiply the values of 
PMN J0948+0022 by a factor ~ 6, which in turn is roughly similar 
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Figure 7. Comparison of the SEDs of the archetypical blazar 1226+023 
(3C 273, grey symbols) and the NLS1 PMN J0948+0022 (red symbols). In 
the case of high-energy 7 rays, the minimum and maximum flux detected 
by Fermi are shown. 



to the ratio of the masses of the central black holes (1.5 x 1O 8 M 
for PMN J0948+0022 and 8 x 10 8 M Q for 3C 273). By shifting up- 
ward the SED of PMN J0948+0022 by this value, the differences 
in the radio-to-X-ray emission become negligible, while those at 7 
rays are emphasized. At these energies, the differences can be ex- 
plained with the different Doppler factor as a consequence of the 
different viewing angle (smaller in the case of this NLS1). 

The basic findings of the present work can be summarized as 
follows: 

• The source reached an observed luminosity of ~ 10 48 erg s _1 
at 7 rays (0.1 — 100 GeV), with day-scale variability and some 
"harder when brighter" spectral behavior. This is the first time that 
such a power is measured from a NLS1. It confirms, beyond any 
reasonable doubt, that NLSls can host relativistic jets as power- 
ful as those in blazars and radio galaxies, despite the relatively low 
mass (1.5 x 1O 8 M0) and the rich environment due to a high accre- 
tion rate (50% the Eddington rate). 

• The coordinated variability observed at all the available fre- 
quencies confirms the typical behavior of relativistic jets, as already 
shown during the MW Campaign 2009. Particularly, the radio emis- 
sion peaked about two months after the outburst at 7 rays, although 
the radio spectral index was already inverted a couple of weeks be- 
fore the peak in the GeV band. The optical flux changed about one 
magnitude (R filter) during the months before the outburst. 

• The morphology at 15 GHz still shows a very compact source 
(~ 10 pc size), despite of the great power released at 7 rays. Com- 
parison with previous observations in 2009 shows that the EVPA 
changed ~ 90° at some time between 2009 July and December 
and this new value is maintained in 2010. The linear polarization 
fraction exceeds 3%, a value greater than that measured in 2009 
« 1%). 

• The comparison with the SED of a typical blazar with a strong 
accretion disk (3C 273) shows that the Compton dominance is more 
extreme in the NLS1. The disagreement of the two SEDs can be 
accounted by the differences in mass of the central black hole and 
Doppler factor of the two jets. 
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